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Abstract 

The paper is devoted to optimization of resonances associated with 1-D wave equations in 
inhomogeneous media. The medium's structure is represented by a nonnegative function B. 
The problem is to design for a given a 6 M a medium that generates a resonance on the line 
a + iM with a minimal possible modulus of the imaginary part. We consider an admissible family 
of mediums that arises in a problem of optimal design for photonic crystals. This admissible 
family is defined by the constraints < b\ < B(x) < 62 with certain constants b\ t 2- The paper 
gives an accurate definition of optimal structures that ensures their existence. We prove that 
optimal structures are piecewise constant functions taking only two extreme possible values b\ 
and &2- This result explains an effect recently observed in numerical experiments. Then we show 
that intervals of constancy of an optimal structure are tied to the phase of the corresponding 
resonant mode and write this connection as a nonlinear eigenvalue problem. 
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1 Introduction 

Recently the increasing interest in loss mechanisms of structured optical and mechanical systems 
has given rise to spectral optimization problems for dissipative models involving wave equations in 
inhomogeneous media, see e.g. (3j [22] and references therein. The question is how to design an 
inhomogeneous medium with very low radiative loss in a given frequency range. The radiative loss 
of energy is closely connected to imaginary parts of eigenvalues of the corresponding non-self-adjoint 
operator, see e.g. [7]. In the paper these eigenvalues are called quasi- (normal) eigenvalues. Naively, 
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the closer quasi-eigenvalues to the real axis R, the less the radiative loss. In the recent numerical 
simulations [T5| IT2] motivated by optimal design problems for photonic crystals, the medium was 
modified by iterative methods with the purpose to move a particular quasi-eigenvalue closer to R. 

An analytic background for spectral optimization problems involving non-self-adjoint operators is 
not well developed. One of the features that make non-self-adjoint spectral optimization problems so 
different from self-adjoint ones is appearance of eigenvalues with algebraic multiplicity greater than 
geometric multiplicity. This leads to a much more complex perturbation theory for eigenvalues. 

The goal of the present paper is to study quasi-eigenvalue optimization problems analytically for 
a 1-D model of a photonic crystal with dissipation at one end. To achieve this aim, we give a rigorous 
treatment of multiple eigenvalues and their perturbations. 

More precisely, the paper is concerned with the eigenvalue problem 

y"(x) + K 2 B(x)y(x) = 0, < x < 1, (1.1) 
2/(0) = 0, (1.2) 
y{l)-xy'{l)/K = 0, (1.3) 

where k is an eigen-parameter and the function B satisfies 

£eL°°(0,l), Q<b l <B(x)<b 2 <oo for a. a. x £ (0, 1). (1.4) 

The spectral problem (jl.ip - (jl.3p comes from the Fourier method applied to the scalar 1-D wave 
equation 

B(x)d?u(x,t) -d%u(x,t) = 0, 0<x<l, t > 0, (1.5) 

equipped with the boundary conditions 

d x u(0,t)=0, d x u{l,t)+d t u{l,t) = 0. (1.6) 

The condition d x u(l,t) +d t u(l,t) = leads to the ^-dependent boundary condition fll.3p and corre- 
sponds to radiative loss of energy to the surrounding medium through the endpoint x = 1. Problems 
(ll.ip - fjl.3p and fll.5p - fll.6p arise in a number of applications. Let us mention spectral problems for 
a 1-D photonic crystal (see e.g. [131 equation (2.26)]) and for an inhomogeneous string damped at 
one end (see e.g. (U [HJJ [2T]). In the first case the physical meaning of the function B is the relative 
permittivity, in the second case B is the density of the string. Note also that b% is always positive 
for optical models, but in the theory of a string, B is allowed to be on a set of positive measure 
[UJ [21]. For studies concerned with the spectral problem f ll.ip -f fL~3]) and the closely related Regge 
problem we refer to [U [191 [ZD El CLU [26] and references therein. 

Eigen-parameters k £ C such that fll.ip -f FOj) has a nonzero solution will be called quasi- 
eigenvalues. The corresponding eigenfunctions are called (quasi-normal) modes. Several other names 
for k are used, sometimes in slightly different settings: dissipation frequencies [HI [21] , resonances 
[T5l IT2] , quasi-normal levels (in the Physics literature) . 

The set of quasi-eigenvalues is denoted by K(B). Quasi-eigenvalues k correspond to monochro- 
matic solutions e lKt if(x, k) of the problem (jl.5p ~ fll.6p . The real part a = Ken of the quasi-eigenvalue 
is the frequency of the monochromatic solution, the imaginary part /3 = Im k is always positive and 
characterizes the rate of decay. 
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The following properties of quasi-eigenvalues are important for the present paper: K(B) is a 
subset of C + symmetric with respect to iR, quasi-eigenvalues are isolated, oo is their only possible 
accumulation point, see e.g. [THJ [2TJ CTJ- 

Let us explain the spectral optimization problem for fll.ip - fll.3p . We take the abstract point of 
view that the problem fll.5p - fll.6p is a mathematical model for a certain device with a structure B(x). 
Assumption f ll.4p defines the family Ad of admissible structures. 

We assume that the device is operated in a particular frequency range [atx, aa], —00 < a± < «2 < 
+00, and denote by Ad[ aij0l2 ] the set of all structures B G Ad such that there exists at least one 
quasi-eigenvalue in this frequency range, i.e., such that K(B) fl {z G C+ : Re 2 G [ai, 02]} 7^ 0- 

The optimization problem under investigation is 

to find B min e Ad [aita2] such that J[ au a 2 ](B min ) < J[ aua2 ](B) for all B G Ad[ aua2 ], (1.7) 
where 

i7[ai,a 2 ](-^) := inf{ImK : Re k G [«i, 0^2] and k G K(B)}. (1.8) 

It seems that the systematic study of eigenvalue's maximization and minimization problems 
associated with self-adjoint elliptic operators was initiated by M.G. Krein [TTJ . In the results of 
[TT] concerned with 1-D self-adjoint problems, extremizers are extreme points of admissible families. 
Krein also proved the same effect for one 2-D optimization problems and conjectured for another 
[TTl Sec. 4. 4]. While there exists an extensive literature on spectral optimization associated with self- 
adjoint elliptic operators (see [HIE] and references therein), there are a very few analytically accurate 
papers on non-self-adjoint spectral optimization problems similar to fll.7p . A possible explanation 
for this fact is that, for self-adjoint problems, eigenvalues move on the real line and do not have 
root eigenf unctions of higher order. This leads to a relatively simple statement of the optimization 
problem and to a relatively simple perturbation theory. Quasi-eigenvalues' behavior is much more 
complex. 

The problems of maximization of the decay rate and of the spectral abscissa for fll.5|) - fll.6p was 
considered in [7J, where existence of the optimal design was proved for a certain class of admissible 
structures from the space I4g' 2 [0, 1] and several estimates on quasi-eigenvalues were obtained. 

In mathematical modeling for photonic crystals, the relative permittivity B is usually discontin- 
uous. That is why the admissible family fll.4p is a reasonable choice. This admissible family with 
hi = 1 was used recently in [T5], where a gradient ascent iterative procedure for optimization of 
the quality factor Q(n) = for individual quasi-eigenvalues k has been developed. Numerical 

computations of [15] were done for the 1-D and 2-D scalar wave equations. Another numerical paper 
[12] is concentrated on the 1-D case, but include an additional coefficient o{x) into the equation 
d x crd x y + K 2 By = and considers various admissible families and discretization techniques. Simula- 
tions of [12] were performed for the case B = 1. It was noticed that for a problem with constraints 
1 = o"i < a(x) < cr 2 = 3 the optimization procedure stopped on a structure a taking only the extreme 
possible values 0^2 . In 1-D simulations of |l5j . the coefficient B also tends to be a piecewise constant 
function taking values 61,2 (see [151 Fig- 1-2 and page 423]). In the author's opinion, figures 3 and 4 
of [15] suggest the same effect in the 2-D case. 

The main results of the present paper are collected in Section [2] (except Theorem 15. 2\ which 
requires more preliminaries). We adjust the existence of minimizer proof of the self-adjoint case [TT] 
to prove that the set of all possible quasi-eigenvalues K(Ad) := {J BeAd K(B) is closed. This easily 
implies that the minimum of j7[ aijQ2 ] is archived and is positive whenever the domain of definition of 
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i7jai,aa] * s nonempty (see Corollary 12. 2p . Then the problem (11.71) can be reduced to the study of the 
case when ct\ = ct 2 = ct. We introduce the function 1(a), a G K, as the minimum of the functional 
J7[ ai a]- It is natural to call the complex points on its graph {(a,Z(a))} ag K = {a + iX(a)} agK optimal 
quasi- eigenvalues, and to call structures B corresponding to k = a + iX(a) optimal structures (for 
the frequency a). (The author does not know whether it is possible that there are non-equivalent 
optimal structures corresponding to a certain frequency q/0). Optimal modes for the frequency a 
are eigenfunctions of (ll.ip -( fl~3l) with corresponding optimal k. 

Theorems 12.41 and 12.51 state that optimal structures are piecewise constant functions taking only 
values b\ and b 2 . For a = 0, we find Z(0) and the corresponding optimal structure explicitly. The 
effect behind Theorem 12.61 is that the intervals where an optimal structure B takes the values b\ or 
b 2 are connected with the arg-function of the corresponding optimal mode. This connection can be 
written in a form of a nonlinear eigenvalue problem. That is, if k,q is an optimal quasi-eigenvalue, 
then the equation 

y" + 4v [ b i + ( & 2 - h) Xc+ (y 2 )] = a.e. on (0, 1), (1.9) 

has a non-trivial solution y Q satisfying boundary conditions (ll.2p - (jl.3p . Here x c+ ( z ) '■= 1 for z G C + , 
and Xr (2) := for 2 G C \ C+. This solution y is an optimal mode corresponding to k . The 
optimal structure associated with k and y is B(x) = bi + (b 2 — &i)x r (yn(x)). 

The results of Section [3] on perturbations of quasi-eigenvalues are preparative for the proofs of 
Theorems 12.41 and 12.51 Perturbations for quasi-eigenvalues of the Schrodinger operator — A + were 
studied, e.g., in [TOj [T], for an abstract approach and more references see [2]. However the proofs 
of our results require more delicate information. Namely, we need analyticity of a quasi-eigenvalue 
as a functional of B and we perform a detailed (though non-complete) study of these functionals 
in vicinity of their singular points. Roughly speaking, these singularities correspond to multiple 
quasi-eigenvalues. The proofs of the main results are based on Lemma 13. 6[ which is essentially 
concerned with two parameter perturbations of a multiple quasi-eigenvalue. Note that the gradient 
algorithm meets obvious difficulties when it encounters a multiple eigenvalue, see the discussion in 
[TB*| p. 425]. An accurate treatment of a multiple eigenvalue requires understanding of its splitting 
picture. Proposition 13.51 and Lemma 13.61 provide a part of this collision and splitting picture. 

Section H] contains the proofs of Theorems 12.41 and 12.51 The proofs are based on the perturbation 
results and on the detailed study of a special solution ip of (II. ip singled out by (p(0) = 1, <p'(0) = 0. 
In Section O we prove Theorem 12.61 and study the interplay of optimal structures B(x) and phases 
&rg(p(x) of associated optimal modes. For the non-degenerate case when b\ > and a ^ 0, this 
interplay is written in an especially transparent form as Theorem 15.21 there exists u G [— it, tt) such 
that B changes its value from bi to b 2 exactly when ip 2 crosses the ray e 1U} M. + and from b 2 to bi exactly 
when ip 2 crosses the ray e lw lR_. 

Notation. <C± = {z G C : ±lmz > 0}, l± = {iGl:±i> 0}, D e (C) := {z G C : \z - C| < e}, 
T = {z G C : |2| = 1}. For £1,^2 sucn that < £ 2 — £1 < tt, Sec(£i,£ 2 ) defines the sector (without 
zero) 

Sec(£i,£ 2 ) := {( G C\0 : argC = £ (mod 2tt) for certain £ G [&,&]}• 

Xe is an indicator function of the set E, i.e., Xe(%) = 1 when x G E, and Xe{%) — when x G" E. 

Open balls in a normed space V are denoted by B e (t> ) := {v G V : \\v — v \\v < e}. For fl C V 
(including the case V = C), v G V, and z G C, let zfl + v := {zv + v : v G Q}. For a function / 
defined on Q C V, /(^) is the image of Q. 
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^c(E)(0> 1) are the Lebesgue spaces of functions with values in C (resp., M); 

W^(0,l):={yeL p c (0,l) : ^£^(0,1), 1 < j < k} 

are Sobolev spaces with standard norms. The space of continuous complex- valued functions with the 
uniform norm is denoted by C[0, 1]. 

d x y, d z ip(x,z;B), etc. denote (ordinary or partial) derivatives with respect to (w.r.t.) x, z, 
etc.; [<9b<£>(x, z; B)](B&) = lim^o <p(x,z,b+c,Bz.)-<p{x,z,b) ^ g ^ e directional derivative of the functional 
tp(x, z\ ■) along the direction B A G L%?(0, 1) \ {0} at the point B G L^(0, 1). 

We write z^ x z^ as n — > oo if the sequences z^/z^ and z^/zf^ are bounded for n large 
enough. 

2 Optimal structures, the definition and main results 

Recall that the set of quasi-eigenvalues corresponding to a structure B (in short, quasi-eigenvalues 
of B) is denoted by K(B). It occurs that K(B) is the set of zeroes of the entire function 

F{z) = F{z; B) := <p(l, z) - i^(l, z) / z, z G C, 

where (p(x, z) = <p(x, z; B) is the solution of the initial value problem 

<%y(x, z) = -z 2 B(x) y(x, z), y(0, z) = 1, 9^(0, z) = 0. 

It is obvious that all modes y corresponding to k G K(B) are equal to ip up to a multiplication 
by a constant. So the geometric multiplicity of any quasi-eigenvalue equals 1. In the following, the 
multiplicity of a quasi-eigenvalue means its algebraic multiplicity. 

Definition 2.1. The multiplicity of a quasi-eigenvalue is its multiplicity as a zero of the entire 
function F(-). A quasi-eigenvalue is called simple if its multiplicity is 1. 

This is classical M.V. Keldysh's definition of multiplicity for eigenvalue problems with an eigen- 
parameter in boundary conditions, see e.g. [23], Sec. 1.2.2-3] and [T9l |2"T] . 

Each quasi-eigenvalue has finite multiplicity. The set of quasi-eigenvalues K(B) is always sym- 
metric w.r.t. the imaginary axis iR, moreover, the multiplicities of symmetric quasi-eigenvalues are 
the same. Note that K(B) may be empty, this is the case for B = and B = 1. These and other 
basic facts can be found in [191 [21] (see also [7]). 

We consider the quasi-eigenvalue problem ( 1 1 . 1 f) - ( T173T) over the following family of structures 

Ad := {B G L£(0, 1) : b x < B(x) < b 2 a.e.}, < h < b 2 < oo, b 2 > 0. 
Proposition 2.1. The set K(Ad) := {J BeAd K(B) is closed and is a subset ofC + . 

The first statement is proved in the next subsection, see Lemma I2T91 The statement K(Ad) C C + 
follows from the well-known fact that K(B) C C + for any B G Ad, see e.g. [HJ l2Tj . 

Now we pass to an immediate corollary, which shows that if K(Ad) has at least one k with 
frequency in the range [01,0^2] (a\,a 2 G R, azi < a 2 ), then a minimizer B m[n for the optimization 
problem fll.7p exists. 

Recall that the functional J[ ai , a2 \ is defined by ( 11. 8ft on the family 

Ad[ aij(X2 ] := {B G Ad : Re k G [ai,a 2 ] for certain k G K(B)}. 
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Corollary 2.2. Suppose Ad[ ai>a2 ] 7^ and consider the problem ( [1. 7| ). Then: 

(i) The functional J[ ai ,a 2 ] achieves its minimum T^ ai OL2 y.= j7[ Qlja2 ](-B min ) over Ad[ aijCe2 ]. 

(ii) The minimum I[ ai _ a2 ] is positive. 

A simple way to check the condition Ad\ aij(X2 i 7^ is to consider 'constant' structures B = b £ 
[61,62] (here and below B = b means B(x) = b for a.a. x £ [0, 1]). For them the quasi-eigenvalues 
are well known. 

Proposition 2.3 (see e.g. [7j). Let B = b be a constant function with 6 > 0. Let {«„} = K(B) be 
the set of corresponding quasi-eigenvalues (taking multiplicities into account). Then: 

(1) 7/6 = orb= 1, then K(B) = 0. 



(ii) Ifb& {0, 1} ; then n n = i^log 



Vb+l 

Vb-i 



+ ^/ n ' %fh>l 77GZ 

+ VE\ n + 1/2, ifb<l ' 



So, excluding the extreme case 61 = 6 2 = 1, we can always ensure Ad[ ai ,a 2 ] taking the 
frequency range [0:1,0:2] wide enough. 

Consider the case 01 = 02 = and introduce the function: 

X(a) = { + °° if M[aM =0 (2 1) 

K J ' \ min{Im/€ : Re k = a and k e K (Ad)} if Ad^j 7^ v ' y 

By Corollary 12. 2[ X(o) > for all o £ M and the minimum in (12 .ip is achieved whenever Adu a i 7^ 
0. Obviously, the minimal value 2J[ Qll a 2 ] °f the functional t 7'[a 1 ,a 2 ](-) is given by min Q:6 [ , 1)Q . 2 i Z(o). 
So problem (11.71) can be reduced to the study of the function X and the properties of structures 
corresponding to quasi-eigenvalues of the form k = o + iX(o). It is natural to call such k and B 
optimal (they are optimal at least for a particular frequency o). 

Definition 2.2. Let X(o) < 00 for certain o £ R. Then: 

(i) k = a + iX(o) is called an optimal quasi- eigenvalue for the frequency o, 

(ii) a structure B £ Ad is called optimal for the frequency o if o + iX(o) £ K(B). 

One can check that the set Ext of extreme points of Ad is 

Ext = {B £ Ad : B(x) £ {61, 6 2 } for a.a. x £ [0, 1]} . 

Denote 

£,(£):= {a; £[0,1] : = &,■}, J = 1,2. (2.2) 

Recall that a function B on [0, 1] is called a piecewise constant function if there exists a partition 
= Xq < X\ < x 2 < • • • < %n < %n+i = 1 such that B is constant on each interval (xj-i,Xj). By 
Ext step we denote the family of piecewise constant functions that belong to Ext (more precisely, the 
family of corresponding classes of equivalence), i.e., 

Ext step = {B £ Ext : Ei(B), E 2 (B) are unions of a finite number of intervals 

after possible correction on sets of zero measure}. 
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Theorem 2.4. Assume that B is an optimal structure for a frequency a G K. T/ien 5 G Ext step , 
i.e., B is a piecewise constant function taking only values b% and b 2 (after possible correction on a 
set of measure zero). 

The proof is given in Section HJ see Corollary I4.10[ Proposition I4.14[ and the proof of Theorem 

For a = 0, we find the optimal quasi-eigenvalue and the corresponding structure explicitly. 
Theorem 2.5. (i) Ifb 2 < I, then 1(0) = +00. 

(ii) If b 2 > I, then 1(0) = log , and the only structure in Ad having the quasi-eigenvalue 
at iX(0) is B = b 2 . 

The proof is given in Section 14.41 We would like to note that statement (i) is equivalent to the 
fact that b 2 < 1 implies K(Ad) PlilR = 0. Under the additional assumption < b\ < b 2 < 1, this fact 
was obtained in Theorem 4.2 (i)] (our proof is completely different). 

In the general case a G R, the intervals where an optimal structure B takes the values &i or b 2 
are connected with the arg-function of the corresponding mode. This connection can be written in 
a concise way as a nonlinear eigenvalue problem. 

Put x c ( z ) '■= 1 when Imz > 0, and x Cj _ (z) '■= when Imz < 0. Consider the nonlinear equation 



y"(x) + K 2 y(x) + (62 - h) Xc+ (y 2 (x))J = a.e. on (0, 1). (2.3) 

Theorem 2.6. Let <b\ <b 2 . Let k be an optimal quasi-eigenvalue for a frequency a6l, Then 
there exists a nonzero solution y G TV^'°°[0, 1] of the nonlinear boundary value problem 112. 3\) . U.2\) . 
U.3\) . Moreover, B(x) = 61 + (b 2 — 6i)x r (y 2 (x)) is an optimal structure for the frequency a. 

Theorem 12.61 is proved in Section [51 where the connection of the rotation of <p 2 (x, k; B) around 
with intervals of constancy of the optimal structure B is explained in details. Note that the solution 
y of Theorem 12.61 is also a mode of the original linear problem (11. ip . (II. 2p . (11.31) . So y(-) = ap(-, k; B) 
with some constant c. In Section [51 we find appropriate constants c = e ld . Generally, 9^0 (mod 27r) 
and ip may be not a solution of (12.31) . ( 11 .211 . (11.31) . However, in some cases y can be taken equal to tp. 

Remark 2.1. (1) Some estimates on 1(a) from above can be easily obtained from Proposition 12.31 
For a certain range of a estimates on 1(a) from below can be obtained using a solution of the direct 
spectral problem for strings of the Krein-Nudelman class |16j . 

(2) Theorem 12.61 does not state that B is a unique optimal structure for the frequency a. The 
author does not know whether it is possible that there are non-equivalent optimal structures corre- 
sponding to certain a / 0. For quasi-eigenvalue optimization problems in classes of Krein strings 
with total mass and statical moment constraints optimal structures are unique for their a 



2.1 Proof of Proposition 12.11 



Lemma 2.7 (Integral form of (ll.ll) - (ll.3p ). A number k G C belongs to K(B) if and only if there 
exists y(x) G C[0, 1] such that 

y(x) = 1-K 2 / (x - s) B(s) y(s) ds, < x < 1, (2.4) 
Jo 

y(l) + ik i B(s)y(s)ds = . (2.5) 
Jo 
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If such y exists, then y(x) = ip(x, K] B). 

Proof. Equality (12.41) holds exactly when y(-) = ip(-,K;B). Using (12. 4p . one can derive (12.51) from 
(jl.3p . And vise versa, equalities (I2.4p - (l2.5p imply that k ^ 0. Indeed, if k — 0, one has y(l) = 0, 
which contradicts (12. 4p . Finally, for k ^ 0, we can rewrite (12. 5p as (II. 3p using (II. ip . □ 

Lemma 2.8. The mapping (z, B) <p(-, z; B) is bounded from C x Lg 3 to Wl'°°[0, 1]. 

Proof. If the families {z^} and {B u } are bounded subsets of C and Lg D (0,l), resp., all functions 
<p(-, z u ; B u ) satisfy the differential inequality < Ci\y(x)\. In turn, this implies (e.g., via the 

Gronwall-Bellman inequality applied to \y\) the statement of the lemma. □ 

Lemma 2.9 (cf. p2] for the self-adjoint case). Assume that there exist sequences {ftn}i° C C and 
{B n }f C Ad such that n n G K(B n ) and K n — > k G C . Then there exists B G Ad such that 
k G K{B). 

Proof. By the sequential Banach-Alaoglu theorem, there exist B G Ad and a subsequence {B n .} 
such that w*— \im.B ri] = B (in weak* topology of L°°). Since K nj — > k, Lemma [2.81 yields that the 
sequence {<£>(•, ft nj ; -Snj)} is bounded in W^'°°[0, 1]. The embedding W^.' °°[0, 1] <e C[0, 1] is compact, 
so there exists a subsequence {m^} of {rij} such that {</?("> ^m 3 \ B mj )} converges strongly in C[0, 1] to 
certain y> G C[0, 1]. This allows one to pass to limits in ( 12 A\i - (12. 51) and to complete the proof using 
Lemma 12.71 . □ 



3 Perturbations and derivatives of quasi-eigenvalues. 

For the sake of convenience, some formal changes in the settings should be done. We extend the 
introduced notation to structures B G L^(0, 1)- For complex- valued B, the statement of the quasi- 
eigenvalue problem and the definition of multiplicities of quasi-eigenvalues remain without changes. 
By ip(x, z; B) and i/j(x, z; B) we denote the solutions of y"{x) = —z 2 B(x)y(x) satisfying 

<p{0,z;B) = d x ,l>(Q,z;B) = l, 0^(0, z; B) = ^(0, z; B) = 0. (3.1) 

Recall that F(z;B) := ip(l,z;B) — id x <p(l, z; B)/z, and that the set K r (B) of quasi-eigenvalues 
of multiplicity r is defined as the set of r-fold zeroes of F(-,B). 

Basic definitions concerning analytic maps in Banach spaces may be found, e.g., in [24J. 

Lemma 3.1. The map (z, B) i— > ip(-, z; B) is analytic from C x L^(0, 1) to W^'°°[0, 1]. Its Maclaurin 
series is 

ip{x, z;B) = l- ip^x; B)z 2 + ip 2 (x; B)z 4 - <p 3 (x; B)z 6 + ..., (3.2) 

<p (x;B) = l, ipj(x;B)= / (x - s) ^_i(s; B)B(s)ds, j G N. 

Jo 

Proof. The W^ft, l]-valued series (13.21) for the solution ip is well-known (see e.g. [H] or [21]). 
It follows from the estimates in [T4"l Sec. 2] (see also [H Exercises 5.4.2-3]) that the series converge 
uniformly on every bounded set of C x L™. So (z,B) i— > ip(-,z;B) is an analytic map on C x 
(see e.g. [Ml Theorem A.2]). □ 
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Lemma 3.2. (i) F(z; B) is analytic on C x L^(0, 1). 
(ii) At quasi- eigenvalues k G K(B), the derivative of F w.r.t. z is given by 



^^ B l (3.3) 



d z F(K] B) = 2 [-/#(1, «; 5) + idMl, k; B)] f ^\s, «; 5)5(s)ds + 

Jo K 

and iae directional derivatives [8 b F(k, B)] {B A ) w.r.t. 5 m the direction B A G L£?(0,1) fry 

[dflFfo S)](S A ) = « [-Kip(l, «; 5) + k; 5)] / <^ 2 (s, k; 5) E A (s) ds. (3.4) 

Jo 

Proof, (i) follows from Lemma 13.11 

(ii) Differentiation w.r.t. B. For any zsC, the functions <^(x, z; 5) and d x (p(x, z\ B) satisfy 



(p(x,z;B) = 1 — z 2 (x — s) B(s) <p(s, z; B) ds, (3.5) 

Jo 

d x (p(x,z;B) = -z 2 B{s) p{s,z;B) ds. (3.6) 



To find directional derivatives [dBf](B A ) and [dgd x ip\(B^), we differentiate these equalities by defi- 
nition using Lemma 13.11 We get 

PX fX 

[dsfix, z; B)](B A ) = -z 2 / (x - s)B(s)[d B <f(x, z\ B)](B A )ds - z 2 / (x - s)B A (s)ip(x, z; B)ds, 

Jo Jo 

(3.7) 

[d B d x <p{x,z;B)]{B A ) = -z 2 [ B(s)[d B ip(s,z;B)](B A )ds-z 2 [ B A {s)<p{s, z; B)ds. (3.8) 

Jo Jo 

It follows from ( 13 .7p that [d B ip(x, z\ B)](B A ) is the solution y(x) of the initial value problem 

y"{x) + z 2 B(x)y(x) = f(x), y(0) = 0, y'(0) = (3.9) 

with f(x) = —z 2 B A (x)ip(x,z;B). Further, (13. 8 p can be rewritten as [d B d x <p(x, z; B)](B A ) = y'(x) 
(as a by-product, we get [d B d x <p](B A ) = d x [d B (p](B A )). Solving (13. 9p by variation of parameters, one 
can find y,y', and, in turn, y(x) — iy'(x)/z. For z / 0, y(x) — iy'(x)/z equals 

f(s)[U/;(x,z;B) - id x ip(x,z;B)/z]<p(s,z;B) - [<p(x,z;B) - id x (p(x,z;B)/z]i/;(s,z;B)\ds. 

Substituting x = 1, f, and z = k G K(B) (so that ip(l, k; B) — \d x ip{l, k;B)/k — 0), we get (13 ,4p . 

Differentiating A3.5\) , 113.6}) w.r.t. z, we see that d z (p is given by the solution y of (13.91) with 
f(x) = —2zB(x){p(x, z\ B) and that d z d x <p = d x d z (p. Hence, for z = k G K(B), 

d z F(n;B) = d z <p(l,K;B) 1 = 2/(1) 1 5 

K K z K K 

id x (p(l, K] B) 



2[-Kij;(l,K;B) + id x i)(l,K;B)} [ V 2 (s, re; B)B(s)ds 

Jo 



K 2 



Using <f(l, re; B) — id x (p(l, re; B)/k = to modify the last term, we get (13. 3p . □ 
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Since the solutions </?(•, z; B) and if)(-, z\ B) are linearly independent, 

at most one of the numbers F(z; B) and [—zi/)(l, z; B) + id x ip(l, z; B)] can be 0. (3.10) 
In particular, 8bF(k, B) ^ if k G K(B). 

Proposition 3.3 (cf. [12] and the discussion in Sec. 5 of [I]). Let k,q G Ki(B q ) (i.e., Ko is a simple 
quasi- eigenvalue) . Then there exist an open ball M p (B ) C -^(O, 1), p > 0, and a unique continuous 
functional k : M p (B ) — > C, such that k(B) G K(B) and k(B ) = k q . Moreover, k is analytic in 
M p (B ) and 

[d B k(B)\(B A ) - — -j— ; — — • (3.11) 

2ac J <£ 2 (s, k o; B)B(s)ds - i(f z {l, k ; B) 

Proof, no G Kx(Bq) if and only if F(ko; Bq) = and 8 z F(ko, B ) ^ 0. The existence of the functional 
Kq(B) with the desired properties follows from Lemma 13.21 (i) and the implicit function theorem 
for analytic maps (see e.g. [24], Appendix B]). To get ( 13. lip , we differentiate F(k(B); B) = 

ip(x) i>(x) 
<p'(x) ip'(x) 

ip(l, k ; B) = id x ip(l,K ;B)/K . □ 



w.r.t. B in the direction B A , and then, use Lemma [3.21 (ii) and the equalities 



1. 



Lemma 3.4. Let P(z, () = z r + hi(()z r ~ 1 + • • • + h r (() be a monic polynomial in z with coefficients 
hj analytic in ( in a neighborhood of ( = 0. Suppose hj(0) — 0, j — 1, . . . , r, and h' r (0) ^ 0. Then 
for ( close to there exist exactly r distinct roots of P(-,() — and these roots are given by an 
r-valued analytic function Z(() that admits a Puiseux series representation 

oo 

Z(() = J2 c j( j/r mth c i = </~ h 'M + (3.12) 

3=1 

(\/C is an arbitrary fixed branch of the multi-function £ 1/V , Cj G C are constants). 

Proof. By the implicit function theorem for multiple zeroes (see e.g. [25| Theorem XII.2]), there exist 
a natural number p < r and a (possibly multi-valued) convergent Puiseux series Z(Q = YlJLi c jC 
such that all its values are roots of P(-, () = for ( small enough. 
Since Z(() = o(l) as ( — > and hj(0) = 0, we see that 

z r (0 = -h r (Q - /v-i(C)z(C) MO^CC) = -CK(o) + o(i)]. 

Therefore Z(() = [-C^(°)] 1/r + o(|C| 1/r ) with h' r {0) ^ for at least one branch of ( 1/r (the equality 
is valid in domains with a cut, e.g., D<s(0) \ e^R + ). This implies p > r, and so, p = r. Further, 
the first term a( 1/r is [-C^(0)] 1/r ^ 0. So Z gives exactly r distinct roots of P(-,C) = 0- Since 
P(-, C) = has at most r roots, we obtain the statement of the lemma. □ 

The following proposition describes splitting of an r-fold quasi-eigenvalue kq under perturbations 
(Ba satisfying [c^F^o; B )](B A ) ^ 0. It is essential that in this case the multiple quasi-eigenvalue 
splits into simple quasi-eigenvalues 'uniformly' in r directions (like roots of z r + C( = 0, C ^ 0). 
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Proposition 3.5. Let k G K r (B ) with r > 2. Assume that B A G Lg 3 satisfies 

I <f 2 (s,K ;B ) B A (s) ds^O. (3.13) 
Jo 

Then there exist open discs 0,5(0), O £ (k ) C C, 5, e > 0, and a convergent in Dj(0) r-valued Puiseux 
series 

fc (0 = «o + E^ ^c 1= ^ ^ ( ^b ) ^° ( 3 - 14 ) 

sitc/i that for any ( G 0,5(0), i/ie r values of k(Q give all the quasi-eigenvalues of B + (B A in D £ (aco) 
and all these r quasi-eigenvalues are distinct and simple. 

Proof. Consider the entire function F(z,Q '■= F(z;B + (B A ) of two complex variables z and (. 
Then k,q is an r-fold zero of the function F(-,0). By the Weierstrass preparation theorem, in a 
certain polydisc B £i (k ) x 0^(0), 

F(Z, C) = [(* - KoY + hiOiz - K y~ l + ■■■ + h r (()]G(z, C), 

where the coefficients hj (the function G(z,Q) are analytic in ©^(/co) (resp., in D £l (« ) x 0^(0)), 
^•(0) = 0, and G(z,() ^ in O £1 (« ) x D 5l (0). 

It follows from flSTTDl . fl3TT3|) . and Lemma O (ii) that [<9 B F(/t ; B )](B a ) ^ 0. Differentiating F 
by (, one gets 

[9bF(«o;S )](5a) = d ( F(K ,0) = ti r {0)G{ Ko ,0). 
On the other side, 9^F(k , 0) = r!G(K ) 0). Hence, 

[gB^ojjojKg^ = rl [dBF(K ;B )](B A ) 

A ) g(«o,o) a;%flo) ^ ' 

So Lemma [33] may be applied to the Weierstrass polynomial (z — k ) t + hi(Q(z — i^o) rl + • • • + h r (C) 
to get the Puiseux series (I3.14p for its zeroes, which are also zeroes of F, and so, are quasi-eigenvalues 
of B + (B A . In this way, we get r distinct quasi-eigenvalues of B + (B A that approach k as ( — > 0. 

Since the obtained r zeroes of F(-,Q are distinct for small (, the standard Rouche's theorem 
argument implies that each of them is of multiplicity 1. □ 

The study of two-parameter perturbations of quasi-eigenvalues requires the following lemma. 
Let us denote 

To := {C=(Ci,C 2 )GM 2 : Ci,C 2 >0 and Ci + C2 < 1}- 

Lemma 3.6. Let Q(z,£i,£z) be a function of three complex variables analytic in a neighborhood of 
the origin = (0,0,0). Assume that is an r-fold zero (1 < r < 00) of the function Q(-,0,0), 
that d(.Q(0) 7^ 0, j = 1,2, and arg9^ 2 Q(0) = arg9^ 1 Q(0) + £ (mod 2n) with £ G (0,7r). Denote 

r\d c .Q(0) 

rjj := — s r Q(o) ' 3 = ^' ^' ^ en f or an V ^ > and £ x in the interval (arg ^/r/i, ^ + arg ^/r/T), t/iere 
exisi a pair (&, C2) G T fl (B*(0) x D 5 (0)) andzeC \ {0} stzc/i toat Ci, C2) = and 
arg z = £1 (mod 27r) . 
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Proof. We give the detailed proof for the more difficult case r > 2. The case when is a simple zero 
of Q(-,0,0) is similar, but much simpler in notation and details. 

It is not an essential restriction to assume that {f- is the analytic in C \ R_ branch of (-) 1//r fixed 
by \/T = 1, and that 

arg«9 Cj Q(0) = (-iyf, arg^-=(-iy|, j = 1,2, (3.15) 

(d r z Q(0) can be placed on M_ by a change of variable in z). 

In this settings, we have to prove that for any £i G (— there exists a sequence {(z n , (^\ (f^)} 

going to and such that arg z n = £i, {C' n ^} C T , and Q(^„, C' n ^) = 0. 

Step 1. An auxiliary triangle and the Weierstrass decomposition. Let us introduce the triangle 

T x = {([1 -6]c,6c) G M 2 : c G [0,1], < 6 1 < 9 < 6 2 < l} , 
with #i, #2 such that 

arg {/^i < < arg {/^, where ^ := (1 - 9j)r]i + 6^772, j = 1, 2. (3.16) 

Clearly, T\ S T , and if a sequence {C' n '} C Ti tends to (0, 0), we have 

Cl n] xC 2 nl x|C N | asw^oo. (3.17) 

By the Weierstrass preparation theorem, in a certain polydisc D £1 (0) x 0^(0) x 0^(0), 

Q(z, C) = P(z, ()R(z, C) with P(z, () = z r + qi {()z r - 1 + ■■■ + q r (C), (3.18) 

where the coefficients qj of the Weierstrass polynomial P (the function R) are analytic in B^OjxB^O) 
(resp., in D ei (0) x D 5 (0) x D 5 (0)), ^(0,0) = 0, and R(z,Ci,G) ^ in D ei (0) x D 5 (0) x B s (0). We 
can suppose that 5 > 1. Indeed, we can always ensure this scaling the variables (1, (2- Now the 
Weierstrass decomposition (I3.18P holds for (Ci, C2) G To- 

Step 2. Asymptotics of zeroes ofQ(z,() for ( G Tx. 

Using the arguments of Proposition 13. 5[ one can show that 

, s d c .Q(0) r\d c .Q(0) 
^(0,0) = ^ = ^^^, , = 1,2. (3.19) 

Assume that there exist z n — > and (Ci^, (2^) — (0, 0) as n — > 00 such that Q(z n ,(^) = and 
{C w } C T x . Then it follows from fl3~T8l) that 

4 = -qr(C [n] ) + g,-i(C w )o(i) + • • • + gi(C N )o(i), n oo. 

Using f)3.19p . (13.171) . and ^(0, 0) = 0, one can show that 

2 \ 1A 

-Ec!%gr(0,0)J [1 + (1)] x|C W | 1/r , n->oo. (3.20) 

.Step ^ Multiplicities of zeroes ofQ(z,() for ( e T t . Let us show that for £ G Ti small enough 
and z small enough, roots of Q(z, () = are simple. 
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Assume the contrary. Then there exist sequences z n and (Ci ; C2 ) as above with the additional 
property that d z Q(z n ,( [n] ) = 0. This implies d z P(z n X [n] ) = 0. Using (KIT} and qj (0,0) = again, 
we see that 

rz'- 1 = q r -i(C [n] ) + ■ ■ ■ + (r - l)ft(C W )<" 2 = K W I 0(1), n -> 00. 

So z n = |C N | 1/(r_1) 0(l). This contradicts fl3T20|) . 

Rescaling d, £ 2 if necessary, we can ensure that all the roots of Q(z, C) = are simple for £ G T\ . 

Step 4- Applying arguments of Proposition 13.51 to the zeroes of the function Q(-, (1 — 0)t,9t) 
with a complex variable r and a fixed parameter 9 G [#1, #2] (for example, = #1), one can produce 
the r- valued Puiseux series 

00 

Z([l - 9}r, Or) = ^ Cj T j ' r with Cl = {/{I - 9) Vl + 9 m ^ 0, (3.21) 
j'=i 

(0j and are from Step 1 and the statement of the lemma). One of the values Zi([l — 9]t,9t) 
of this function can be chosen if we place the branch {/t instead of the multi-function r 1 ^ in the 
Puiseux series (13. 2 1 p . Step 3 and the implicit function theorem for simple zeros imply that Z\ can 
be extended to an analytic on T\ \ {0} function Zx(Ci, (2)- For each 9 G [9i, 62] and sufficiently small 
r G (O,ei(0)], e x (9) > 0, the function Z X {[1 - 9}t,9t) is given by the Puiseux series 1ET2I]) with t 1 ' 1 ' 
replaced by its branch {/t and coefficients Cj depending on 9. (However, we do not know if e{9) is 
uniformly separated from 0, and so, we use other arguments to study the asymptotics of Z x ). 
Comparing ( I3.2ip with ( I3.20p and using (" 13 . 1 7[) , one can see that 

^i(Ci,Ca)= V / Ci^i + C2^ 2 [l + o(l)], C^0, (gTl (3.22) 

This implies that for ( G T\ small enough, Z X (C) 7^ and, due to (13. 15ft . axgZx(() is a continuous 
function with values in (— tt/2, 7r/2). On the other side, ( 13.201) and ( 13 . 161) imply that for r small 
enough, argZ 1 ([l — 9x\t,9xt) < £1 and argZiQl — 9 2 }t, 9 2 t) > £1. So for each r G (0,e 2 ), e 2 > 0, 
there exists 9 G (#1, 9 2 ) such that arg Zx([l — 9)r, Or) = £1. This completes the proof. □ 



4 Proofs of Theorems 12.41 and 12. 5 L 

4.1 Optimal structures are extreme points of Ad; the case a ^ 0, b\ > 0. 

We start from the case when b x > 0. Note that this assumption is satisfied in the optimization 
problem for photonic crystals [151 H2]- F° r the degenerate case bx = 0, some details of the proof 
require modifications (see Subsection 14.31 below) . 

Lemma 4.1. Let B G L\(0, 1), B(x) > a.e., and z 2 g R. Then: 
(i) ip(x, z;B)^0 for all x G [0, 1] . 

(*) 9 ~W^T ? R f° r aU xe 

(Hi) For any £ G (— ir, ir] the set {x G [0, 1] : <p(x, z\ B) G e^R + } is finite. 

Proof, (i) Let <p(xx,z;B) = 0. Then x x > (since ip(0,z;B) = 1), and so, <p(x,z;B) is an eigen- 
function of the self-adjoint boundary value problem — y" = XBy, y'(0) = y(xx) = 0. Hence the 
corresponding eigenvalue A = z 2 is real, a contradiction. 
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(ii) Let h = ^^-b) e M for x i > °- Then ^( x > z \ B ) 

is an eigenfunction of the self-adjoint 
problem — y" = XBy, y'(0) = 0, y'(xi) — hy{x\) = 0. So z 2 G R, a contradiction. 

(iii) Let := {i G [0, 1] : </?(x, z; 5) G e^R + } be infinite. Then it includes a convergent sequence 
of distinct points {x n } nS N C ft. It follows from statement (i) and ip G C[0, 1] that xo := limx n also 
belongs to Q. Then agzJzwio) g e^R, and so, 9 x <^(xo) G e^R. If xo > 0, this contradicts statement 

So Xq = 0. Assume z 2 G C_ (the case z 2 G C + is similar). Then there exists a neighborhood of 
(— z 2 ) lying in a sector separated from R, i.e., there exist an interval [£i,£2] C (0,7r) and £ > such 
that 

B £ (-z 2 ) C{(GC\0 : argCG [6,6]}=: Secfe^). 

By continuity of (p(x) = <p(x,z;B), there exists 5 > such that {— z 2 <£>(x) : a; G [0,5]} C D £ (— z 2 ). 
Since B > a.e., we have —z 2 B(x)(p(x) G Sec(£i,£ 2 ) f° r a - a - x e [0,5]. Using (I3.5p . one can show 
that <f(x) G 1 + Sec(£i,£ 2 ) f° r ah x G (0,5]- The intersection of the ray e^R + with 1 + Sec(£i,£ 2 ) 
is either empty or separated from the point 1 (note that ^ Sec(£i,£ 2 )). Using <£>(0) = 1 and 
(p G C[0, 1], we can choose 5 so small that <p(x) ^ e^M.+ for all x G (0, 5]. So xo = is not a limit 
point of Q, a contradiction. 

□ 

Lemma 4.2. Let B G Ad ; 2 G C. Then the set of directional derivatives [dsF(z; B )](Ad — B ) := 
{[d B F(z; B )}(B A ) : B A + B Q e Ad} is convex. 

The lemma follows from the fact that Ad — B := {B A : B A + B Q G Ad} is convex. 

Lemma 4.3. Let Bq G Ad \ Ext, k,q G K(Bq), and Re fco 7^ 0. JTien an interior point of the set 
[d B F(K ;B )}(Ad-B ). 

Proof for the case b\ > 0. Since -Bo £ Ad \ Ext, there exist E\ > such that the set Q := {x G (0, 1) : 
61 + £1 < £ (x) < 62 — £1} is of positive Lebesgue measure. Using Lemma 1431 one can show that 
there exist two sectors Sec(£i,£ 2 ) an d Sec (£3, £4) such that the sets 

ill ■= {% e Q : ip 2 (x,K Q ;B ) G Sec(£i,£ 2 )} and fi 2 := {x G Q : y? 2 (x, k ; -Bo) e Sec(£ 3 ,£ 4 )} 

are of positive measure, and 

Sec(£ 1 ,£ 2 ) n Sec(£ 3 ,£ 4 ) = 0, Sec(£ 1 ,£ 2 ) n (-l)Sec(£ 3 , £4) = 0. 

This, k G K(B ), and (13 . 1 j) imply that the complex points (1 and £ 2 defined by 

Cj ■= k [-k iI)(1,k ;B ) + id x ip(l,Ko;B )] / y? 2 (s, k ; B )ds, j = 1,2, 

are linearly independent as vectors in R 2 . Due to Lemma 13721 (ii). Q = [8bF(ko, -B )](xo.) . 

Taking B A j(x) = exuAx), e G (0,£i), j = 1,2, we see that B A j G Ad — £> (by the definition of 
the set fl). So the complex intervals Ci( — £ i> £ i)j J = 1> 2, are subsets of [9bF(k ; -Bo)](Ad — 5 ). 

Since £1 and £2 are linearly independent, is an interior point of the convex hull of the intervals 
£i(— £1, £1) and £ 2 (— £1, £1). Lemma I4~2l concludes the proof. □ 



15 



Proposition 4.4. Let B Q G Ad and Kq G K(B ). If zs an interior point of the set 
[8bF{hq\ B )](Ad — Bq), then there exists 0i > such that k — G X(Acf). 

Proof. The quasi-eigenvalue Ko of So has a finite multiplicity r (see e.g. j2TJ Ej)- Recall that 
k G K r (B ) implies d r z F(K ; B ) ^ 0. 

Since is an interior point of [9b.F(/cq; Bo)](Ad — B ), there exist Ba,i,Ba,2 G Ac? — £> and a 
branch of multi- function £ 1 / r such that 



By Lemma 13.21 (i), Q(z, Ci, C2) : = E(k + z; B Q + Ci-Sa,i + C2-5a,2) is an entire function of three 
complex variables if, 0, and £2- Obviously, /to is an r-fold zero of Q(-, 0, 0), and due to the properties 
(i)-(iv), we can apply Lemma T3.6I to the function Q. Lemma T3.6I implies that there exist (Ci, C2) G To 
and z G C \ {0} such that Q(z, Ci> C2) = and argz = — n/2 (mod 27r). So 2f = — /3i > 0, and 
k - 0i G K(Bo + Ci^a.i + C2£ A , 2 ). 

Since Ad — Bq is convex and contains 0, the structure Bq + Ci-E?a.i + C2-Ba,2 belongs to Ad. □ 

We see that Lemma [4.31 and Proposition 14.41 imply the desired statement (for the case 61 > 0). 
Corollary 4.5. Let a ^ 0, B G Ad, and a + iZ(a) G #(-B). T/ien 5 G Ext. 

4.2 Piecewise constancy of optimal structures for a 7^ 0, 61 > 0. 
Definition 4.1. Let 61 7^ 6 2 and B G Ext. 

(i) We say that x G [0, 1] is a switch point of B if for any e > both the sets Ej(B) n (a; — e, x + e), 
j = 1,2, are of positive measure (see (12.21) for the definition of Ej(B)). 

(ii) We say that a switch point Xq is singular from the left (right) if for any e > both the sets 
Ej(B) n (x — e,x ) (resp., Ej(B) fl (x ,x + e)), j = 1,2, are of positive measure. Otherwise, x is 
called regular from the left (resp., right). 

(iii) A switch point Xq is said to be regular if it is regular both from the left and from the right. 

It is easy to notice the following. 

Lemma 4.6. Let B G Ext have a switch point Xq singular from the left (right). Then there exists a 
sequence {x n }^° of (distinct) switch points converging to Xq from the left (resp., right). 

This implies easily the next statement. 

Lemma 4.7. Let B G Ext. Then the following statements are equivalent: 

(i) B G Ext s tep. 

(ii) B has a finite number of switch points. 



(i) The complex numbers 7]j : 



r! [d B F(K ;B )](B AJ ) 



, j = 1, 2, are nonzero. 



dr z F(K ;B ) 
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(Hi) All switch points of B are regular. 

Let us define a quasi-argument function qarg : C \ {0} — > T by qargz = e iarg2: . 
Note that for z G C + \ iR and B G Lr(0, 1), 5 > a.e., Lemma [4.11 implies that, for a suitable 
branch of the argument function, 

arg<yj 2 (x, z\ B) is continuously different iable on [0, 1] and d x arg<£> 2 (x, z; B) ^ for x G (0, 1]. 

(4.1) 

In the next lemma, we keep such a choice of the arg-function. 

Proposition 4.8. Let < b±. Let Bq G Ext have a switch point xq G (0, 1). Let kq G K{Bq) \ iR ; 

£(x) = arg</> 2 (x, n ;B ), and & := aTg(K [-K ip(l, k ; B ) + id x i)(l, k ; B )}) 
(£i is well defined due to A3.10\) ). Then the set of quasi- arguments of directional derivatives 

A(k ; So) := qarg [d B F(K ; B )} ([Ad - B ] \ {0}) (4.2) 

contains at least one of the semi-circles {e 1 ^ ^ 1 ^ : s G (0,7r)}, {e 1 ^ ) 4 ^ 14 " 5 ! : s G (— 7r,0)}. 

More precisely, assume that the measures of the sets (x — e, x ) fl Ej^B) and (x , x + e) n E^, 
where ji,j 2 G {1, 2} an<i j\ ^ j 2 , are nonzero for all e > 0. Then 

(i) in the cases {^'(^o) < 0, ji = \, j 2 = 2} and { (,'( x o) > 0, ji = 2, j 2 — 1 }, the set A(k, ; Bq) 
contains { e 1 ^ x ^ + ^ +s ^ : s G (0,tt)} ; 

(^mJ in the cases {£'( x o) < 0, j\ = 2, j 2 = 1} and {£'(x ) > ; ji = 1, j 2 = 2}, the set A(k, ; B ) 
contains {e^ ^ 14 ^ : s G (-7r,0)}. 

Proof. We consider the case when £'( x o) < and ji = 1, j 2 = 2. Arguments for all other possible 
cases are similar. Then there exist two sequences {L n }^°, {i? n }^° of subsets of (0, 1) such that 

(LR1) All L n and R n are of positive measure. 
(LR2) L n c (x - i x Q ), C (x ,x + J), n G N. 
(LR3) -Bo(^) = bi for x G L„, B (x) = b 2 for x G 

The assumptions on L n and i? n imply that {b 2 — b\)\L n and (&i — b 2 )xR n are in Ac? — B . So 
the set [d B F(K ;B )](Ad - B ) is convex, contains r? L ,n := [d B F(K ; B )\ ([b 2 - &i]xO, ^7-R,n := 
[9b-F(/«o; -Bo)] ([&i — ^2]Xfi„)) and, obviously, contains the point 0. 

Now finding from (13. 4p the arguments of f]L,n, VR,n f° r large n, one can prove the statement of the 
lemma. Indeed, for any e > 0, taking n large enough, we can ensure that £(x) G (£(^o), £(^o) + e) 
for x G L n , and that £(x) G (£(x ) — e, £(x )) for x E R n (the assumption £'( x o) < is used here). 
From this and fl3.4[) we see that argr] i n = £1 + (,L,n (mod 27r), where 



:= arg 



(62 - 61) / y2 2 (s, K o; B ) ds 



e (£(*o),f(zo) + <0 



Similarly, taking into account b\ — b 2 < 0, 

arg^ >n = £1 + 6j„ ( mod 2vr), ^ n G (£(x ) + vr - e, £(x ) + vr). 

So the quasi-arguments of the convex hull of 0, i]L,ni and r] R:n cover at least the set {e 1 ^ ) 4 ^ 14 ^] : 
s G [e, 7r — e]}. Moving e — )■ 0, we complete the proof. □ 
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Proposition 4.9. Let B Q G Ext /iai>e a singular switch point x G [0,1]. £e£ k G K(B ) and 
Re k 7^ 0. T/ien is an interior point of the set [8bF(ko; B )](Ad — B ). 

Proof for the case b\ > 0. By Lemma |4.6[ there exists a sequence of distinct switch points {xj}^ 
converging to x . Let = &rg(p 2 (x, k ; B ) as in Proposition 14.81 By Lemma H~T1 (iii), there 
exist j and n such that e^ x ^ 7^ e^ Xn \ e^ x ^ 7^ — e l ^ Xn \ and Xj,x n G (0, 1). Applying Proposition 
14.81 to the switch points Xj and x n , we see that A(k ; B ) contains two distinct semi-circles. Since 
[dsF^Ko] B )](Ad — B Q ) is convex and contains 0, this implies the proposition. □ 

Corollary 4.10. Let a 7^ ; B G Ad, and k = a + iX(a) G K(B ). Then B G Ext s t ep . 

Proof for the case bi > 0. By Corollary 14. 5[ B G Ext. Assume that B g" Ext step . Then, by Lemma 
14.71 i?o has a singular switch point. Propositions 14.91 and 14.41 imply that there exists K\ G K(Ad) 
such that Re^! = a, but Imux < T(a). This contradicts the definition of T(a). □ 

4.3 The case 61 = 0. 

In this subsection we complete the proof of Corollary 14.101 considering the case b\ = 0. 

While the most part of the proof remains the same, some changes have to be done in Lemma 14.11 
and Proposition 14.81 The case b\ = is degenerate in the sense that the set E\(B) (where B equals 
0) may have a positive measure and include an interval [0, xq) with xq > 0. 

Define 

a 1 := sup{x G [0, 1] : B = a.e. on [0, ai]}. (4.3) 
So B(x) = a.e. on [0, ai], and (if a\ < 1) 

B is not equivalent to on any interval (ai,x ) with xq>cl\. (4.4) 

Statements (ii) and (iii) of Lemma [4.11 are not valid on [0, ai]. Lemma [4.11 can be adjusted in the 
following way. 

Lemma 4.11. Let B G L^(0, 1), B(x) > a.e., and ai < 1. Let z 2 g" E. Then: 
(1) (p(x, z; B) ^ for all x G [0, 1] . 

(iii) For any £ G (— 7r, 7r] the set {x G [ai, 1] : f{x, z; B) G e^R + } is finite. 

Proof, (i) Obviously, <p(x, z; B) — 1 for x G [0, a x ]. Let <p(xi, z; B) — for sx > a x . Then <yj(x, 2; 5) 
is an eigenfunction of the nontrivial self-adjoint problem defined on the interval [ai,xi] by (11. ip and 
the boundary conditions y'( a i) = v( x i) = ( see e -g- [HE]). The corresponding eigenvalue A = z 2 
is real, a contradiction. In the same way we get statement (ii). 

(iii) We put Q := {x G [di, 1] : <p(x,z;B) G e^R+} and assume that Q is infinite. Then, as 
before in Lemma 14.11 Q has a limit point xq G Q, and combining this with statement (ii) of the 
lemma, one can show that Xq — a±. 

Since B = on [0,ai], <p satisfies ip(x,z;B) = 1 — z 2 f ai (x — s)B(s)y(s)ds (in particular, 
<f(ai, z; B) = 1). This and f|4.4p implies that x 7^ a\ in the same way as before. □ 
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Proof of Lemma \4-3\ in the case b\ = 0. The proof requires only the change of Lemma H~T1 to Lemma 
14.111 Indeed, one can see that Q := {x G (0, 1) : b\ + E\ < B (x) < b 2 — £1} is a subset of (eti, 1]. So 
Lemma [4.111 can be applied to show the existence of the sets Qi and Q 2 . D 

Now the proof of Corollary 14.51 is complete. 

The following changes are needed in connection with Proposition 14.81 As before, for a suitable 
branch of the argument function, argcp 2 (x, z; B ) is differentiable. However, d x arg (p 2 (x, z; B ) ^ 
only for i£ (ai, 1] (this follows from Lemma [4.111 (ii)). 

Proposition 4.12. Let B G Ext have a switch point xq G (ai,l). Then all the assertions of 



Proposition 4-8 hold true 



Since x G (a 1; 1), the proof of Proposition 14.81 works without changes. Note that, for B G Ext, 
the definition of a\ implies that Bq has no switch points in [0, ai). 

After all these modifications, the proofs of Proposition 14791 and Corollary 14. 101 requires no 
changes for the case b\ = except the use of Proposition 14.121 Now the proof of Theorem 12.41 
for a 7^ is complete. 

4.4 The case a = 0. 

In this subsection we will study quasi-eigenvalues on the ray ilR + . In this case, the analysis is simpler 
since 

<f(x, z] B) and ip(x, z\ B) are real when z G iM+ and B G L^(0, 1). (4-5) 

Lemma 4.13. Let Bq G Ad and kq G ilR+ fl K(Bq). If Bq ^ b\ and Bq ^ b 2 , then there exist 
S A ,+,5 A ,- G Ad- B such that [d B F(K , B )](B A:+ ) > and [8 b F(kq, B )](B A -) < 0. 

Proof. It follows from fl3TT0l) . Kty . and k = i(3, > 0, that 

K [-KoV(l, «ol S ) + id x i>(l, k ; B )] G K \ {0}. (4.6) 

The assumptions of the lemma imply that there exist sets Qi,Q 2 C [0, 1] of positive measure and 
e > such that B (x) > b\ + e for x G Qi and B (x) < b 2 — e for x G Q 2 . Put B A j = (— l^e^n., 
j = 1,2. Then, B A j G Ad — B , j = 1, 2. Since the real continuous function ip(-, k ; Bq) has at most 
finite number of zeroes on [0, 1], we see that 

i ,i 

(p 2 (s, k ; £0) B At i(s) ds < 0, / ip 2 (s, k ; B ) B At2 (s) ds > 0. 

Jo 

Combining this with (13. 4p and (14. 6p . one gets the statement of the lemma. □ 

Proposition 4.14. Let Bq G Ad and kq G iR+ fl K(Bq). If Bq ^ b\ and Bq ^ b 2 , then there exists 
0i > such that k - i/5 x G K(Ad). 

Proof. We give a proof for the case when kq is a quasi-eigenvalue of multiplicity r > 2. The arguments 
for the case of a simple quasi-eigenvalue are simpler in details. 

So d r z F(n ; B ) ^ 0. It follows from gS} that F(i0; B ) G E when G R+. Hence, d r z F(K ; B ) G 
i r M \ {0}. By Lemma [4.13[ we can choose B A G Ad — Bq such that 

/ r\[d B F(K ;B )}(B A ) \ 
arg ^ dr F{Ko . Bo) J = -™/2 2tt). 
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So one can choose yf- such that 



<■' r\[d B F(K Q ;B )](B A ) 
arg ,/ dr F{Ko . Bo) = (mod 2n). (4.7) 



Applying Proposition 13.51 we see that for £ > small enough, the r branches of the Puiseux series 
(I3.14p give all the quasi-eigenvalues of B + (B& that tend to k as ( — > 0. Let us choose in (13. 14ft 
the branch of (-) 1//r analytic on M + and such that l l ' r = 1. Then it follows from (14. 7p that there 
exist a quasi-eigenvalue ki(() of B + (B&, C > 0, with asymptotics &i(C) = K o — ic|C| 1//r + °( I C I 1//t ") , 
c > 0. The other branches of k(() have the same asymptotics for ( > with constants c ^ R + . Since 
K(B + C-Ba) is symmetric w.r.t. iR, &i(C) stays on iIR + for ( > small enough. 

Summarizing, we see that the structures B + (B& belong to Ad for ( G [0, 1], and that for ( > 
small enough, one of quasi-eigenvalues of B + (B& may be written in the form k — i/3(C) with 
/5(C) > 0. ~ □ 

Proof of Theorem \2.5l Let k,q = iX(0) G iM + be a quasi-eigenvalue of B e A<i. Then Proposition 
14.141 implies that either B = b\ or B = b 2 . Proposition 12.31 completes the proof. □ 



5 The proof of Theorem 12.61 : nonlinear eigenvalues and re- 
strictions on switch points. 

Let B be an optimal structure for a frequency a^0 and let n = a + il(a) 6 K(B) be a corresponding 
optimal quasi-eigenvalue. By Theorem 12.41 B is a piecewise constant function taking only the values 
bi and b 2 - In other words, B has at most finite number of switch points {xj}™ =1 where B changes its 
value either from b\ to b 2 or inversely from b 2 to b\. Note that the values of B at the switch points 
and the endpoints x = and x = 1 are not important for the quasi-eigenvalue problem. 

We will use the notation of Proposition 14.81 with the continuous in x branch £(x) of the multi- 
function &rg(p 2 (x, k; B) fixed by £(0) = 0. 

Lemma 5.1. If an optimal structure B is constant on an interval (xi,x 2 ), then \£,(xi) — ^,{x 2 )\ < n. 

Proof. Assume the contrary, i.e., — > 7T- Then (13. 4p implies that the set of quasi- 

arguments A(n; B) defined by (14.21) contains an arc of length greater than ir. So is an interior point 
the convex set [ObF^ B)](Ad—B). Proposition l4.4l implies that k is not an optimal quasi-eigenvalue, 
a contradiction. □ 

Consider first the non-degenerate case < b\ < b 2 . Since k G C + \ iR, the function £ is 
continuously differentiable on [0, 1] and £'( x ) 7^ for x G (0, 1]. So £,'(x) keeps its sign on (0, 1]. In 
fact, for x G (0, 1], 

£'(£)< if Re k > 0, and g(x) > if Re n < 0. (5.1) 
(see (14. ip and the proof of statement (iii) of Lemma [4.11 for details). 

Theorem 5.2. Let < 61 < b 2 . Let B be an optimal structure for a frequency a G R\ {0} and let k 
be a corresponding optimal quasi-eigenvalue. Then there exists u> G [— tt, 7t) such that, on the interval 
(0, 1), B(x) changes its value from bi to b 2 exactly when (p 2 (x, k; B) intersects the ray e lw lR + and B 
changes its value from b 2 to bi exactly when <£> 2 (x, k; B) intersects the ray e lul M._. 
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Proof. We assume that < on (0, 1] (arguments for the case > are similar). 

Case 1. Assume, first, that B is constant on (0, 1). Then Lemma |5 . 1 1 implies that £(1) > —it and 
the statement of the theorem is valid both with u = —tc and u = 0. 

Case 2. Assume that B has only one switch point x\. Then, in the case when B changes its 
value from b\ to b 2 at x±, Lemma T5 . 1 1 yields the statement of the theorem with u = £(sci), and in the 
opposite case with uj = + it. 

Case 3. Consider the case when there are at least two switch points. Let the set {xj}™ =1 of switch 
points be naturally ordered < x\ < X2 < ■ ■ ■ < x n -\ < x n < 1. To be specific, assume that B 
changes its value from b\ to b 2 at the first switch point x\. 

By Lemma I5TTI £(xi) G [— 7r, 0). Assign uj = £(xi). From Proposition 14.81 and the assumption 
< 0, we see that the set A(k; B) contains the semi-circles {e 1 ^^'^ 1 " 1 "^ : s G (0, 7r)} for odd j 
and the semi-circles {e 1 ^^"^ 1 ^ : s G (— 7r,0)} for even j. Since k is an optimal quasi-eigenvalue, 
all these semi-circles coincide (see the proofs of Propositions 14.91 and Corollary 14.101) . This means 
that 



Thus, B may change its value from b\ to b 2 only at the points where £(x) —00 (mod 27r) and from 
b 2 to 61 only where £(x) = u + it (mod 2ir). On the other side, Lemma [5. II implies that B indeed 
changes its value each time when £(x) = u (mod 2%) or £(x) = u + ir (mod 2ir). This completes the 
proof. 

In the case when B changes its value from 6 2 to b\ at x\, analogous arguments produce u = 



The proof of Theorem \2.6\ in the case b\ > 0. Let a = 0. Then an optimal quasi-eigenvalue k exists 
only if b 2 > 1. The optimal structure is B = b 2 , and the corresponding mode ip(-,K; b 2 ) is a real- 
valued function with a finite number of zeroes (see Section \AA\) . Hence, for arbitrary 9 G (0,tt/2), 
the function y = e l9 (p is a solution of the nonlinear eigenvalue problem (12.31) . (II. 2p . (II. 3p . 

Now consider a^0 and take uj as in Theorem 15.21 Then, to obtain a solution of the nonlinear 
problem (pO]) . (fl~2|) . (Q|) . one can put y(-) = e i(7r ~ w)/ V(-, «; 5) in the case Re/t > 0, and y(-) = 
e _lw / 2 <^(-, K] B) in the case Re n < 0. This fact follows immediately from (15. ip and Theorem 15.21 □ 

If 61 = 0, some technical complications arise since ( 15. ip is valid only on (ai, 1], while i\x) = 
for x G [0, ai]. Recall that [0, a{\ is the greatest interval of the form [0,a;o] such that B(x) = a.e. 
on [0, Xq], see Section H~3l for details. 

The proof of Theorem \2.6\ in the degenerate case bi = 0. In the cases when a = or ai = 0, the 
proof is the same as in the case bi > 0. 

Consider the remaining possibility when a^0 and B(x) = on (0,ai) with < a\ < 1. Then 
ai is the first switch point of B and (p(x, k; B) — 1 for x G [0, ai]. With no loss of generality, assume 
additionally that (,'( x ) < on (a 1; 1] (that is, we assume Ren > 0). 

Let us show that, on the interval (ai, 1), the optimal structure B changes its value from to b 2 
exactly when (p 2 (x, k ; B) intersects M + and from b 2 to exactly when <p 2 (x, Kq] B) intersects R_. 

Indeed, under the above assumptions, the arguments of Proposition 14.81 applied to the first switch 
point X\ = ai imply that A(k; B) contains {e'K 1+s l : s G (0,7r)} (note that £(ai) = 0). Since (15. ip 



{fa) 



C(xi) (mod27r) for odd j, 
£(xi) + 7T (mod 2ir) for even j. 



(5.2) 
(5.3) 



CM + 7T. 



□ 
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holds true on (ai, 1], we can deal with all the other switch points x%, £3, ... in the same way as in 
the proof of Theorem 15.21 As a result, we obtain that the semi-circle {e t ^ 1+s l : s G (0, 7r)} generated 
by the first switch point and all the semi-circles 

{ e iK(*j)+ei+«] : s e (o ;7r )} for odd j > 1, {e i[ ^ )Hl+s] : s e (-n,0)} for even j, (5.4) 

generated by the other switch points coincide. This yields the desired statement. 

Now one can take y(x) = e~ 17T / 2 Lp(x, k; B) and check that y is a solution of the nonlinear problem 
O, (H2D, (IT31). 

Note that in the case £'(2) > on (ai, 1] (i.e., when Res < 0), these arguments produce y(x) = 

if(x,K]B). □ 
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